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We investigate the flux penetration patterns and matching fields of a long cylindrical wire of 
circular cross section in the presence of an external magnetic field. For this study we write the 
London theory for a long cylinder both for the mixed and Meissner state, with boundary conditions 
appropriate for this geometry. Using the Monte Carlo Simulated Annealing method, the free energy 
of the mixed state is minimized with respect to the vortex position and we obtain the ground state 
of the vortex lattice for N = 3 up to 18 vortices. The free energy of the Meissner and mixed state 
O ' provides expressions for the matching fields. We find that, as in the case of samples of different 

, geometry, the finite size effect provokes a delay on the vortex penetration and a vortex accumulation 

Oh' in the center of the sample. The vortex patterns obtained are in good agreement with experimental 

results. 
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PACS numbers: 74.60.Ec, 74.60.Ge 



I. INTRODUCTION 



^3 

O, 

Modern technology has made possible to fabricate superconducting samples of small size, like films of thickness 
less than the London penetration length as well as superconducting wires of radii of order of this length. This has 
I ' arisen the interest to study again the geometrical or size effects in superconductors. The size problem was already 
^ , studied a long time agotra but recently it has been reconsidered and some studies of transport curmnts, magnetization, 
magnetic moment, reversibility lines, and flux penetration in finite size samples have been madeoil3. In particular, the 
penetration of magnetic flux in finite size superconductors has attracted the attention of physicists in the last years 
• due to the fact that in finite samples the geometrical effects can produce important modifications in the critical state 
and in the resulting vortex distribution. As an example we can cite the results of Zeldov et al.Lj in thin superconducting 
strips which have been found a delay in the penetration of the vortex lines and a vortex accumulation at the center 
of the sample due to the geometrical barrier effect, but other studies at this respect has been madecl. In this way, our 
main interest in this work is to study specifically the size effects in the vortex lattice and in the matching fields (the 
minimum field for a new vortex penetration) for an infinite superconducting wire of circular cross section. 

We use the London theory to study the formation of vortex lattice in the mixed state of a long superconducting 
cylinder. London theory is valid in the limit of low induction (fields well bellow the upper critical field H C 2), where 
most of the experiments can be performed to observe vortex lattice. This theory fails in the limit of small length 
scale. This breakdown of the theory is originated in the fact that the finite size of the vortex cores is neglected. 
Thus, London theory is not suitable to treat the self-energy of a single vortex line. However, what really matters to 
determine the shape of the ground state of the vortex lattice is the interaction energy between vortices on different 
sites. For an infinite sample, vortices interact one with each other via a two-body potential. Nevertheless, if the 
surface is taken into account, additional terms which describe the interaction of the vortex with the surface come 
in. In this work we develop the theory to describe such interactions in a long cylinder. In addition, we use London 
theory to determine the free energy of the Meissner state of a long superconducting wire. For the determination of the 
ground state of the vortex lattice we use the Monte Carlo Simulated Annealing minimization method. For the first 
time, we obtain the ground state lattice pattern starting from an arbitrary configuration. This procedure is different 
than the usually used where the free energy of the vortex lattice for several predetermined configurations is calculated 
and pick the lowest one. 

The paper is outlined as follows. In Section [n] we determine the magnetic field of an arbitrary distribution of 
vortices in the mixed state and its corresponding induction (spatial average of the local magnetic field). In addition, 



we find both the London and Gibbs free energy. In Section III we repeat this calculation to the Meissner state. In 
Section IV and |v| we analyze the matching fields and the vortex lattice patterns. 
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II. MIXED STATE 



In what follows we develop a theory for the mixed state of a long superconducting wire. Our starting point is 
the London equation. This equation is obtained from the second Ginzburg-Landau equation by assuming that the 
superconducting order parameter is a constant throughout the whole space, that is, it neglects variations of the order 
parameter inside the vortex-core. The London approximation is valid provided that the Ginzburg-Landau parameter 
k = A/£ 3> 1; here £ is the coherence length, and A the London penetration length. In addition, the vortices, whose 
size is of the order of A , may overlap but not the vortex-cores. In cylindrical polar coordinates (r, 4>) the London 
equation for the local magnetic field h = hz is given by 

x2 fd 2 h ldh 1 d 2 h\ , _ ^ . 

where $o is the quantum flux, r^, is the position of the z-th vortex inside the cylinder, and <5(r) = 5(x)5(y) is the two 
dimensional delta function. Here we are assuming that the vortices are straight lines. Therefore, a 3D problem is 
reduced to a 2D one. We will solve this equation subject to the following boundary conditions 

h(a,(t)) = H 
dh 



., ; =0, (2) 

where a is the radius of the cylinder. The first condition assures that outside the sample the field is uniform and is 
precisely the external field H ; the second one states that the perpendicular component of the current vanishes at the 
boundary r = a, that is, the Cooper pairs cannot jump out of the sample. 

To solve (0) we use the Green's function method. Assuming for the Green function the boundary conditions 
G(a,(/>,r' ,()>') — 0, G(r,(j),r' ,<j)') continuous at r = r', dG(r,(f),r' ,<f>')/dr discontinuous at r = r' and assuming that 
both G and h are periodic in <p, one obtains 



h(r', <t>') = $oV G(n, fa,r', <f>') - HaX 2 f 

Jo 



dG 
dr 



(3) 



To proceed we need to solve the equation for the Green's function. The method we use to find this function is 
outlined in Ref. [ |ll| except by the fact that there the Green's function is associated with the Poisson equation and 
boundary conditions are taken at infinity. One has 

G(r, 0, r', 41) = ^ [KoQt - r'\/X) - a(r, 0, r' , </>')] , (4) 



where 



a{r,far'J)= £ cos[m(0 - 4')} I m (r/X)I m (r' /X) , (5) 
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where I m and K m are the modified Bessel functions. 

We are now in a position to find the local magnetic field. Substituting (^) into (^) and using the identity 
I m (x)K' m (x) — I' m {x)K m {x) — — - (where the primes on I and K stands for the first derivative with respect to 
x) to develop the second term for the field expression one obtains 

Hr, 0) = ^E [^d r - r «l/ A ) - ^ & r *>k)] + H J^j ■ ^ 

The London (Hclmholtz in the thermodynamic context) free energy contains basically two contributions. One is 
the energy stored in the field and the other one is the kinetic energy of the supercurrents. The London free energy 
per unit length is 



F _ 1 
L ~ 8^ 



pa /-27T 






/ / drd 
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dhy X 2 (dh 
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x h(r i ,cf> i ) + — — d<p [ — ) , (7) 
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where on going from the first to the second line we have used the London equation (|]J) and the boundary conditions 
(||) and the periodicity of the field. Here L is length of the system. 
Now, the London free energy can be evaluated by introducing (^jj) into (Q). This yields, 



F = f <Pp V 
L \4ir\J 



^[ifodri-rjl/A)-^^ 



H\ ahiflJX) 
~2 A 7 (a/A) 



(8) 



where H = J ff/($ /47rA 2 ). 

To obtain the equilibrium configuration of the vortex lattice the Helmholtz free energy is not convenient because 
the calculations involve a fixed magnetic field H . Therefore, it is necessary to perform a Legendre transformation to 
obtain the Gibbs free energy. The Gibbs free energy (in units of volume) is given by, 



BH 

47T 



(9) 



where T — F/AL, with A = tvo 2 , and B is the induction which is the spatial average of the local magnetic field, 

f 



B = -| J d 2 rh. 

The evaluation of this integral is tedious but straightforward. We obtain 

N<f> Q 2na\H h{a/\) $ 1 
~ ~ 4~~ + A Jo (a/ A) ~ ~AI (a/X) *f 

Finally we obtain for the Gibbs free energy 



(10) 
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(12) 



Notice that the first term is the vortex self-energy and the second one describes the repulsive interaction between the 
vortices (the bulk term) . The third term describes the attractive interaction between the vortices and image vortices 
located outside the sample. The effect of this interaction is to push the vortices close to the surface. The argument of 
the fourth term represent the flux $/$o = (1 — ^o( r )/^o( a ))- The first term on this argument represents the repulsive 
interaction between a vortex and the magnetic field that penetrates the sample surface and it push the vortices to 
the center of the sample. The second term in this argument represent the vortex magnetic energy. The fifth term is 
the Meissner state energy (see next Section). The competition between the vortex- vortex image interaction and the 
interaction between the vortex and the surface field represents an energy barrier that the vortex has to overcome to 
be able to enter the sample. When the external field is bellow the matching value, the vortex-vortex image interaction 
is more important and a new vortex is not able to enter. When the field is increased up to or above the critical value, 
the vortex can overcome the surface energy barrier and enters into the sample. On the other hand, we want to point 
out that the bulk interaction is invariant under any translation but this symmetry is no longer valid for the whole free 
energy because of the presence of the various interactions. However the system is still invariant under any rotation 
because it depends only on the angle difference. 



III. MEISSNER STATE 



In the Meissner state, although we have no penetration of vortex lines we have penetration of magnetic field near 
the surface. In a semi-infinite superconductor, for instance, the external field penetrates exponentially over a distance 
A. For a superconductor with cylindrical geometry, the field inside the sample is given by the London equation, 

_A 2 (— + -— \ +h-0 (13) 
\ dr 2 r dr J 
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The boundary condition is 

h(a) = H . 

The solution for this equation with the appropriate boundary conditions is 

-WA) 



h{r) = H 



I (a/X) 



Using eq. (0) the London free energy per unit length can be written as 

F = / $o \ 2 (&V a h{a/\) 
L \4tt\J \ 2 I AJ (a/A) ' 



(14) 



(15) 



(16) 



which is the second term of (||). Notice that if a 3> A, near the surface the field has an exponential behavior 
h ~ H(a/r)e {r - a V x like in a semi- infinite superconductor. 
Using equation (10) we find for the induction 
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2iza\H h(a/X) 
A J (o/A) ' 



(17) 



which is the second term of (|l l|) . 

The Gibbs free energy can be calculated by substituting (|l(]) and ( |l7|) into @. One has 



Q = - 



$0 

4ttA 2 



.4 



a Ii(a/ A) 
A I {a/ A) 



(18) 



which is the forth term of ( p^ ) . The previous results comprises the framework for the discussion of the superconducting 
properties of a long cylinder. In the next Sections we examine the matching fields and the vortex lattice patterns. 



IV. MATCHING FIELDS 



The lower critical field H c \ , defined as the lowest external field enough to have penetration of at least one vortex 
line, is determined assuming that at the phase transition from the Meissner state to the mixed state, the Gibbs free 
energy per unit length, i.e., QA, has the same value. In what follows, S stands for the mixed state and M for the 
Meissner state. One has, 



M 



ABmH c i 

47T 



AT, 



AB s H cl 

47T 



(19) 



where (AFm, ABm) are given by ([l6 17) and (AFs, AB$) by (js|,|Tl"[) respectively. Introducing these equations into 
( |l9| ) we find for the lower critical field 



4ttA 2 



In k 



K (a/X) 
/o(a/A) 



Io(a/\) 



(20) 



Here we have taken the center of the cylinder as the equilibrium position of a single vortex line in (|3|) and ( |1 1| ) . In 
this case, only the m = term survives in the sum for er (cf. equation @)). By taking the limit of a — > oo in (|2(i| ) 
we recover the well known result = (<I>o/47rA 2 ) hiK. In Fig. [j] we plot the difference AH = H c \ — H^!f in units 
of (<f>o/47rA 2 ) which shows that the smaller the values of a the larger the value of the lower critical field. In this way, 
this result shows clearly that the size effect, provokes a delay in the first flux penetration as found in experiments on 
other finite systems of different geometrylm 

As a consequence of the energy barrier generated by the finite size of the sample we have a delay not only for the 
first vortex line penetration but for the subsequent lines too. In this way, at low fields (near the lower critical field), 
we have a well defined critical field for each new penetration (matching field). Different from the bulk case, where the 
induction increases continuously with the external field, here the induction increases by steps as showed in Fig. |[ A 
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new vortex enter in the sample only when the energy is enough to overcome the surface energy barrier. As we go to 
higher magnetic fields the induction approaches the bulk case. 

The matching fields for each configuration of vortices (H s n — H s jy /($o/47rA 2 ), N = 1, .., 18; with H sl = H c i) are 
calculated using the same procedure used to obtain H c \, i.e. equating the Gibbs free energy of the configurations 
with N and N + 1 vortices, Gn — Gn+i- In this way we obtain a transcendental equation in H (because the radii 
implicitly depends on H) which can be solved iteratively. However, each time the free energy is calculated for a 
given magnetic field, this free energy must be calculated with the vortices in the equilibrium position. Thus, for each 
iteration the free energy must be minimized. The iterative wjarJk. was performed using the Secant method and the 
minimization using Monte-Carlo Simulated Annealing methoca£3. For a/A = 10 the values obtained for H s n can be 
seen in Table |. 



V. VORTEX LATTICES 



The usual procedure to find the ground state of the vortex lattice is to assume some particular geometry and then 
evaluate the Gibbs free energy. The configuration corresponding to the lowest value is the one supposedly the most 
stable vortex lattice. Other authorsESEj have been used the method of images to determine de vortex [Gopijiguration. 
In the present work we follow a different procedure. Using Monte Carlo Simulated Annealing methodi-JO, we start 
from an initial configuration chosen randomly and we let the vortex lattice evolve towards the global minimum. The 
energy minimization is made using different initial configurations and different seeds for the random number generator. 
In this way we obtain different annealing schedules and assures that the system goes to the global minimum. This 
procedure has been usually avoided by many authors because the computational time strongly increases with the 
number of vortices. However, with the advent of very fast computers this is no longer a major problem. 

As the minimization method may require a large number of evaluations of the Gibbs free energy, we must find 
some effective manner of calculating <j(r, (j>, r' , <j>') (cf. equation ([5])), which involves a sum of infinite terms. We will 

consider a ^> A, but still finite. Within this approximation 7 m (o/A) ~ , 1 e°/ A , and K m (a/X) ~ \ ?wT\ e ~ a ^ X - 

W27ra/X y Za / A 

One has, 

a(r, 0, r', <f>') = ^- 2a / A / (|r + r'\/X) , (21) 
where we have used the following identityB 

^ cos(m(j))I m (x)I m (y) = h(\/x 2 +y 2 + 2xycoscf)) . (22) 

m— — oo 

The Gibbs free energy can be simplified to 

g = 



4>n \ 2 A 2 



4ttA 2 J A 
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V Ma 



hin/X) ( hY a hja/X) 

(a/A) ' I 2 J X I (a/X) 



(23) 



where r; = — r-j. 

To analyze the vortex patterns for each number of vortices we choose arbitrarily the magnetic field in the middle of 
the interval between H s n and H s n + i. The minimization of the free energy was performed again using the Simulated 
Annealing method. As we have observed in Section [o], the system is invariant under any rotation. Therefore, we can 
fix one of the vortices along the x axis so that the minimization procedure will involve 2N — 1 variables. We have 



done this for N = 3 up to N = 18. The vortex patterns for N = 3 up to N = 11 are illustrated in Fig. VI and 
the complete numerical data for N = 3 to 18 is shown in Table Q. As it can be seen from the figure and table, the 
vortices arrange themselves in quite simple geometries. For example, for N — 3 an equilateral triangle; for TV = 4 a 
square; for ]V = 5a pentagon; for N = 6 a pentagon with a vortex at the center; for N = 7 a hexagon with a vortex 
at the center; and for N = 18 a vortex at the center, a first ring of six vortices forming a hexagon and a second ring 
of eleven almost equally spaced vortices. The results show a clear tendency of the inner vortices to form a hexagonal 
lattice, however, the external not. For N — 3 to 9 the radii of all the vortices in the ring coincide, however, for N > 9 
the vortices in the outer ring show small radii fluctuations (no more than 7%). For example, as we can see on Fig. 
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VI for the case of N — 10, we have two vortices in a symmetric position related to the center of the sample and 
an external ring with 8 vortices. The vortices closer to the central ones have bigger radius than those closer to the 
empty space between them. This is easy to understand based on the competition between the different interaction as 
we pointed out at the end of section II. When the external filed is larger than H s iq, the repulsive interaction of the 
vortices with the field that penetrates the surface is more important than the vortex-vortex image interaction. This 
interaction competes with the repulsive interaction between the vortices which is more important between the central 
vortices and the vortices of the outer ring that are closer to them. This result agree with the result of Ternovskii et 
al.tj. They showed for a semi-infinite plane, .that the vortex lattice is distorted near the surface. 

As long as a decade ago Yarmchuk et alx3 determined experimentally the vortex patterns up to 11 vortices in 
superfluids and Campbell and ZiflO made a very extensive numerical study of these patterns. Despite the interaction 
between vortices in a superconductor is screened over on scale A, and in supcrfluid it is logarithmic on all scales, there 
is indeed a clear analogy between them. The patterns found here are in excellent agreement with the experimental 
results found in Ref. [ [l^] . The vortices does not show only similar patterns but show the same tendency to vortex 
accumulation at the center of the sample. The-same tendency has been observed in experiments with finite size 
superconducting samples of different geometriesncl and this suggest that this tendency is proper of finite samples 
independent of the particular geometry. Our results are in good agreement also with the results of Ref. [ ^| , if we 
compare the ring structure but, two differences we can point out. First, as they supposed predetermined ring structures 
their vortex patterns does not show the detailed structure of the outer rings as our more accurate calculations. Second, 
their vortex patterns does not show the tendency of vortex accumulation at the center of the sample. 



VI. SUMMARY 



In summary, by using London theory we have studied the size effects on a long superconducting cylinder both in the 
mixed and Meissner state. We have determined numerically the ground state of the system and obtained the vortex 
patterns, the induction and matching fields for N = 3 up to N = 18 vortices. Our results show a clear tendency of 
vortex accumulation in the center of the sample, delay in the penetration of the flux lines and irregularities in the 
distribution of the vortices closer to the surface, due to finite size effects. Those effects have been found before in 
experiments with samples of different geometries showing, apparently, that this is a consequence of the finite size of 
the samples and not of the particular geometry. 
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FIG. 2. Induction (73) for TV = 1 to 18 vortices as a function of the external magnetic field (iff) in units of ($o/47rA 2 ). 
For each N vortices configuration, the induction is calculated in the [7ff s jv, ff s jv+i] field range. Here we have used a/A = 10, 
A = 200nm and k = 40. 
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Fig. 3/Venegas/Sardella 



FIG. 3. Vortex patterns for different values for N = 3 up to JV = 11. The parameters used here are a/A = 10, A = 200?im, 
and k = 40. The distances between adjacent circles is 2 (in units of A). 
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TABLE I. Here No represents a vortex at the center of the sample, N 1 and N2 are the number of vortices on the first and 
second ring and Ri and R2 are the respective radii of these rings. From 3 to 9 vortices, all the radii coincide up to the last 
digit shown in the table. Above 9 vortices, the internal radii (Ri) coincide but, because of the fluctuations in the value of the 
external radii, we have assumed R2 as the mean value. The seventh column is the energy associated to the terms that depend 
on the vortex positions and the last column represents the values of the critical fields for each configuration of N vortices. The 
critical fields were determined with a precision of 1CT 7 and the energies of 10 10 , respectively, though these two quantities are 
quoted with seven decimals. Here we have used the same parameters as in Fig. H and IVD. 
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